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Abstract 


A design method has been developed by which an airfoil with a substantial amount 
of natural laminar flow can be designed, while maintaining other aerodynamic and geo- 
metric constraints. After obtaining the initial airfoil’s pressure distribution at the design 
lift coefficient using an Euler solver coupled with an integral turbulent boundary layer 
method, the calculations from a laminar boundary layer solver are used by a stability anal- 
ysis code to obtain estimates of the transition location (using N-Factors) for the starting 
airfoil. A new design method then calculates a target pressure distribution that will 
increase the laminar flow toward the desired amount. An airfoil design method is then 
iteratively used to design an airfoil that possesses that target pressure distribution. The 
new airfoil’s boundary layer stability characteristics are determined, and this iterative pro- 
cess continues until an airfoil is designed that meets the laminar flow requirement and as 
many of the other constraints as possible. 
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1.0 Introduction 

Since Orville Wright first flew in December of 1903, there have been considerable 
attempts to find methods to reduce the drag of airplanes. A reduction in drag means that 
airplanes can operate more efficiently by using less fuel, which results in lower operating 
costs and smaller, quieter engines. Also with the reduction in fuel consumption comes the 
ability to produce aircraft with longer ranges and bigger payloads. 

In the 1930’s, it was found that longer runs of laminar flow over an airfoil resulted in a 
lower profile drag and that favorable pressure gradients contributed to prolonged laminar 
boundary layers (ref. 1). Using these ideas, pressure distributions having the pressure 
minimum located near the position of desired transition were sought. Once the desired 
pressure distribution was found, an airfoil with that pressure distribution was then derived, 
using theoretical techniques such as Theodorsen’s method (ref. 2), and tested. The NACA 
1-6 series airfoils are examples of airfoils that were designed in this manner (ref. 3). This 
was the birth of attempts to achieve long runs of natural laminar flow (NLF) to reduce air- 
plane drag. 

In the 1960’s, a new method for creating long runs of laminar flow was utilized. Now 
called laminar flow control (LFC), this method achieved laminar flow through suction 
holes located at selected spanwise stations on the wing. There are two results of boundary 
layer suction. First of all, boundary layer suction thins the boundary layer and lowers the 
effective Reynolds number. Secondly, boundary layer suction changes the boundary layer 
profiles. The changes that result contribute to boundary layer stability, which results in 
longer runs of laminar flow (ref. 4). Since then, this technology has been further 
researched and several airfoils using LFC have been developed (refs. 5,6). 
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Although the benefits of LFC are tremendous, especially in three-dimensional flows, 
the physical application of an LFC system to a wing causes several problems. One prob- 
lem is the increased weight that the system adds to the airplane. Since the aircraft weight 
is increased, a trade-off must be made to get the true benefits of the LFC system. Another 
problem of LFC is contamination from insects and icing. Often times, insect remains or 
ice on the surface trip the boundary layer, reducing the efficiency of the LFC system. 
Another possible problem arises if a mechanical failure occurs and the system on one 
wing does not work properly. In this case, lift will be lost on the wing and its drag will 
increase, causing unwanted rolling and yawing moments. If these problems with LFC 
systems can be overcome, perhaps the true benefits of the technology can be experienced. 

Within the past fifteen years, due to the advances of the modem-day computer, even 
more research has gone into the benefits of and the methods for achieving NLF. Although 
most of this research has been applied to designing airfoils (refs. 7-11), some research has 
been done on designing fuselages for NLF (refs. 12-14). The methods implemented for 
designing airfoils with long runs of NLF seem to be mainly trial and error methods using 
linear stability theory to assess the effect of changing the pressure distribution, although 
there have been some optimization methods developed for axisymmetric airplane ele- 
ments (ref. 12). Thus, the design methods employed for modifying airfoil pressure distri- 
butions often required extensive knowledge and experience. 

In this paper, a constrained design method is presented for modifying an airfoil’s shape 
such that long runs of NLF can be achieved. The design method uses an Euler solver cou- 
pled with an integral turbulent boundary layer method and a laminar boundary layer solver 
to calculate the velocity and temperature profiles at each airfoil station. A boundary layer 
stability analysis code is then used to find the stability of the boundary layer in terms of N- 
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factors, which are the logarithmic amplification of the Tollmien-Schlichting (T-S) waves. 
After calculating a target N-Factor distribution that forces transition to occur at the desired 
location, a new method is used to calculate a target pressure distribution that is closer to 
meeting the NLF constraints, while maintaining several aerodynamic and geometric con- 
straints. Once the new target pressure distribution is calculated, an airfoil design method 
is used to design an airfoil that has that pressure distribution. This process is iterated until 
an airfoil is designed that meets the desired NLF, aerodynamic and geometric constraints. 
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2.0 Overview of the Airfoil Design Method 

The flowchart shown in Figure 1 demonstrates the process by which the final NLF air- 
foil is iteratively designed. The first module on the flowchart uses the Euler solver dis- 
cussed in Section 2.1 and the turbulent boundary layer method discussed in Section 2.2 to 
calculate the pressure distribution of the current airfoil. This pressure distribution is then 
analyzed by the laminar boundary layer solver discussed in Section 2.3 to calculate the 
boundary layer velocity and temperature profiles. The stability analysis code discussed in 
Section 2.4 then uses this data to calculate the N-Factor distribution for the current airfoil 
and pressure distribution. Using the current N-Factors and pressures, the Target Pressure 
Design module, which is discussed in Chapter 3, calculates a target pressure distribution, 
from which an airfoil can be designed using the flow solver and the airfoil design method 
discussed in Section 2.5. However, the target pressures often need to be modified while 
designing the new airfoil. These modifications are made by the Modify Target Pressures 
to Enforce Constraints module, which is discussed in Chapter 4. 

Chapter 5 is included to show some results of the NLF airfoil design method. Airfoils 
were designed for glider, commuter and subsonic transport aircraft, which covers a wide 
range of Mach numbers, Reynolds numbers, and airfoil thicknesses. 

While Chapters 3 and 4 discuss the newly developed components used in this project, 
the current chapter is devoted to discussing the existing codes that have been coupled 
together to design new NLF airfoils. The computer used was a Silicon Graphics Indigo2 
workstation with an R4000 processor; all CPU times mentioned pertain to this computer. 

2.1 Euler Solver 

The Euler solver used is called GAUSS2 (General purpose Approximate factorization 
Upwind Scheme with Shock fitting, 2-Dimensional version), and was written by Peter 
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Hartwich (refs. 15-17). The code solves the two-dimensional, compressible, non-conser- 
vative Euler equations on a structured mesh. 

GAUSS2 uses an upwind method that constructs a finite-difference scheme based on 
the theory of characteristics. Since the upwinding technique used is the non-conservative 
split-coefficient matrix (SCM) method, the code is fast and efficient. The dependent vari- 
ables used are the speed of sound, and the two Cartesian velocity components, all of which 
are required in any upwind scheme. Entropy is chosen as the fourth dependent variable, 
which reduces the energy equation to a simple convection equation. 

Away from shocks, fully one-sided, second-order accurate spatial differences are used 
to update the solution. The shocks are resolved using a floating shock fitting method that 
allows the shock to float between two grid points. Across the shock, the Rankine-Hugo- 
niot relations are explicitly used to update the solution using the speed of sound and 
entropy variables. 

The solution is advanced in time using a time-implicit operator containing block-tridi- 
agonal matrices for the two-dimensional Euler equations. For calculating transient flows, 
the second-order accurate Crank-Nicholson time differencing is used. For the steady-state 
calculations used in this study, the first-order accurate Euler-backward time discretization 
is used due to the quick convergence rate that results from its better damping properties. 

As mentioned previously, GAUSS2 is a fast flow solver. Typically only 500 CPU sec- 
onds are required to calculate the pressure distribution of the initial airfoil at the design lift 
coefficient. 

As shown in Figure 1, GAUSS2 must also be used each time that the airfoil design 
method is used to design a new airfoil. Typically, only 60 iterations through the airfoil 
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design method are needed to design an airfoil that possesses the desired target pressures. 
Approximately 500 CPU seconds are required to complete these 60 iterations. 

Required inputs for GAUSS2 include the current airfoil, the angle of attack, and the 
free-stream Mach number. Outputs of the solver include the pressure distribution on the 
surface of the airfoil, the location of any shocks on either surface, and the wave drag asso- 
ciated with these shocks. 

2.2 Turbulent Boundary Layer Method 

The compressible turbulent boundary layer method used is a modified version of the 

integral method developed by Stratford and Beavers (ref. 18). The following seven equa- 
tions are taken from reference 18. With M being the local Mach number, an equivalent 
flat plate distance is defined as 



This distance represents the length over which a boundary layer growing on a flat plate 
would acquire the same thickness as the real boundary layer at that given location and 
Mach number (ref. 19). 

For free-stream Reynolds numbers Re between 1 and 10 million, the method 
expresses the boundary layer thickness as 

8 = 0.37X O8O R<’'°- 20 , (2) 

the momentum thickness as 

0 = 0.036(l +0.10AT 2 )“° 7 V 80 Re' 0 ' 20 , (3) 

and the displacement thickness as 
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6* = 0.046 ( 1 + 0.80M 2 j°' 44 X 0 '* 0 /?e' 0 ’ 20 (4) 

For free-stream Reynolds numbers between 10 and 100 million, the method uses the fol- 
lowing relations: 


8 = 0.23X 0 - 833 /?,- 0167 


e = 0.022(1 -t-o.ioAf 2 ] -070 ^ 833 ^ 
8* = 0.028^ 1 + 0.80A/ 2 J 044 X°' 833 Re 


- 0.167 


- 0.167 


(5) 

(6) 
(7) 


In order to find expressions for the preceding equations that are a good approximation 
throughout the entire range of Reynolds numbers, the equations above have been modi- 
fied. The expressions used in the current application of the boundary layer method are 


8 = 0.276X°' 82 fo“°' 18 

(8) 

0 = 0.027(l - i-O.IOM 2 ] -0 ' 70 ^' 82 /?^ 0 ' 18 

(9) 

8* = 0.034(l +0.80M 2 J 044 X°‘ 82 Re“ 018 

(10) 


The original method of Stratford and Beavers does not include a calculation of the vis- 
cous drag. As a result, the method of Squire and Young (ref. 20) is used for this purpose 
and is implemented after calculating the characteristics of the turbulent boundary layer as 
described above. This method extrapolates the momentum thickness at the trailing-edge 
of the airfoil to infinity. 

In addition, the original method of Stratford and Beavers does not include a criterion 
for predicting the location of turbulent separation. Therefore, the method of reference 19 
is included for this reason. This method is based on the pressure gradient parameter 


( dC\ i'2( Vo.io 

F = C p{ X ~d?) I 10 XRe ) ’ 


( 11 ) 
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where X was defined in Equation 1. 

The displacement thicknesses calculated using Equation 10 arc used to calculate an 
effective airfoil to account for viscous effects. Withy being the airfoil station, the upper 
surface of the effective airfoil is calculated using the relation 


i - 1 




( 12 ) 


where y e j u is the upper surface of the effective airfoil, yj u is the upper surface of the 


current airfoil, k is a relaxation factor (typically 0.80), and (6 *) y „ and (5*) ; u are the 

upper surface displacement thicknesses of the previous and current iterations, respec- 
tively. To calculate the lower surface of the effective airfoil, the relation implemented is 

y e,j, I = yj. i - K ( 5 *)y, I - ( 1 - K ) ( 5 *) j~l 1 (13) 

It is this effective airfoil that is analyzed by GAUSS2 in order to calculate a pressure dis- 
tribution. 

2.3 Laminar Boundary Layer Solver 

The laminar boundary layer solver used is BL3D, which was written by Venkit Iyer 
(ref. 21). This code is a quick, compressible, three-dimensional solver with fourth-order 
accuracy in the wall-normal direction. It is applicable to attached laminar flows where the 
perfect gas and boundary layer assumptions are valid. 

Since the solver is being used for two-dimensional applications, the x- and y-momen- 
tum, continuity and energy equations are solved step-wise starting at the stagnation point 
of the flow. The method is iterated at each station until the desired convergence is 
obtained. After the solution at each station has converged, the derivative of the velocity in 
the wall-normal direction at the surface is checked to see if the flow has separated. At 
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each station, the velocity and temperature profiles, their first and second derivatives, and 
several boundary layer edge conditions are written out for the stability analysis code. 
These are the only outputs that are needed from this solver. The inputs required are the 
Mach number, Reynolds number, and the airfoil with its pressure distribution. 

2.4 Stability Analysis Code 

The stability analysis code used in this method is the COSAL (Compressible Stability 
Analysis) code written by Mujeeb Malik (ref. 22). In order to obtain the stability proper- 
ties of three-dimensional compressible boundary layers, COSAL solves the eigenvalue 
problem by solving an eighth-order system of differential equations. For a two-dimen- 
sional case, the problem reduces to a sixth-order system of equations. 

Small disturbance theory was used to obtain the basic equations for the linear stability 
analysis of parallel-flow compressible boundary layers. A set of five ordinary differential 
equations result from the compressible Navier-Stokes equations and small disturbance 
theory. These equations are composed of one first-order continuity equation, three sec- 
ond-order momentum equations, and one second-order energy equation. This system can 
be reduced to a set of eight first-order ordinary differential equations. 

Using a finite difference method, COSAL solves the system of basic equations. In 
order to do this, the eigenvalues are initially obtained through a global eigenvalue search 
since there are no guesses available. Then, once a guess is obtained, a quick local eigen- 
value search is used to continue solving the equations. 

COSAL uses temporal stability theory which assumes that the disturbances grow or 
decay only with time. As a result, the frequency o> is assumed complex, 

to = + (14) 


9 



while the wave numbers C, and £ are assumed real. The disturbances are said to grow if 
©j. > 0 , and decay if G) ( < 0 . 

With a complex group velocity defined as 



an N-factor used for transition prediction is calculated using the relation 



where a represents the arc length along a curve on the surface being analyzed. 

COSAL can use four different methods to integrate to find the N-Factor at each x-loca- 
tion. These methods include the envelope method, the fixed wavelength and orientation 
method, the fixed wavelength and frequency method, and the fixed orientation and fre- 
quency method. The method used in the current application of COSAL is the envelope 
method which requires the real frequency co r to be specified, and then maximizes the 

growth rate co ( with respect to the wave numbers C and \ . 

Inputs to COSAL include the real frequency, an initial wave angle and the boundary 
layer profile data from the laminar boundary layer solver. In this application, the most 
important output from COSAL is the N-Factor distribution. 

When obtaining the stability characteristics of an airfoil, usually a minimum of 10 fre- 
quencies must be analyzed by COSAL. For supercritical cases, perhaps as many as 20 or 
25 frequencies must be analyzed. At 200 CPU seconds per frequency, COSAL requires at 
least 2000 CPU seconds, but perhaps as many as 5000 CPU seconds, for each airfoil 
design iteration. This makes the use of COSAL the most expensive aspect of the method. 
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Approximately 70% of the time required to use this NLF airfoil design method is spent in 
COSAL. 

2.4.1 Calculation of Analysis N-Factor Distribution 

As mentioned above, the envelope method in COSAL is used to obtain the N-Factor 
distribution of the particular pressure distribution being analyzed. This method requires 
the real part of the frequency to^ and an initial wave angle y to be specified in order to 

obtain a unique N-Factor distribution by maximizing the growth rate co ( with respect to 
the wave numbers £ and ^ . An infinite number of frequencies a> r occur in nature in the 

flow over an airfoil. Therefore, in theory, in order to calculate the exact N-Factor distribu- 
tion for the pressure distribution being analyzed, COSAL would have to be used to calcu- 
late an N-Factor distribution for every frequency that exists in nature. 

Since this is not realistic, the pressure distribution is analyzed for a range of frequen- 
cies for each airfoil station j from the stagnation point (J = 1 ) to the laminar separation 
point 0 = k ). The lower and upper frequencies of this range, denoted by co r mjn and 

co r max respectively, and ny, the number of frequencies to be analyzed within this range, 
are specified. From n y, the differential frequency Aoo r is calculated as 


co — co 

a r, max r, min 

Aco r = 

n r \ 


(17) 


This means that the initial frequencies that COSAL will analyze are co r min , 


co 


. + Aco , co + 2Aco , ..., co - Aco , and co 

r, min r ’ r, min r' ' r, max r' r,max 

The N-Factor at a particular x-location on the airfoil will differ with frequency. For a 
constant wave angle, the N-Factor at the x-location will increase as the frequency is 
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increased. However, at some critical frequency co r crJf , the N-Factor will be a maximum 

and, as the frequency is increased beyond that, it will start to decrease. This is demon- 
strated in Figure 2 in which a typical N-Factor is shown over a range of frequencies at a 
point on an airfoil at a constant wave angle. In general, co r crit is different for each x- 
location on the airfoil. As a result, even though each of the frequencies in the specified 
range have been analyzed by COSAL, this does not guarantee that u> r crit for each x-loca- 

tion has been determined. 

Therefore, a method had to be developed that would ensure that to r crit for each airfoil 
station was obtained. A flowchart of this method appears as Figure 3. The process begins 
by defining two variables: Nj and co^. With j denoting the airfoil station, Nj is the array 

of maximum N-Factors, while is the array of critical frequencies corresponding to the 
N-Factors of Nj. Initially, both arrays Nj and o Oj are set equal to zero at every j. 

When COSAL is called for the first time, the pressure distribution is analyzed at a fre- 
quency of co r min and the wave angle \\f . The N-Factor of (i) r min at station j (defined as 

N u j) is compared to N } at ;. If N (0j >N j , then AT is replaced with the value of N aj , 
and (Hj is replaced with 0) r min . Then, the same process is repeated for each of the fre- 
quencies in the frequency range, as shown in the flowchart. After this has been done, Nj 
contains the first estimate for the array of maximum N-Factors, and to contains the first 
estimate for the critical frequencies. 
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However, if any frequency in array (Oj is equal to to r min , then the N-Factor at that sta- 
tion j (found in array Nj ) is not really the maximum. Using the data in Figure 2 as an 
example, if the specified frequency range was 30000 Hz to 40000 Hz, then the frequency 
found in cor for this particular station j would be 30000 Hz since the N-Factor at this sta- 
tion is the largest N-Factor in this range. But looking at the figure, the N-Factor at 30000 
Hz is not really the maximum. As a result, frequencies below to r mjn must be analyzed in 

order to reach the maximum N-Factor for station j . This is accomplished by extending 
the frequency range by changing the value of 0) r min to to r min - Aco r . After analyzing 

this frequency, array c Oj is checked to see if it is equal to co r min - Aco r for any j . If it is, 
then a» r min needs to be modified again. If 0) r min - Aco r does not appear in array o D- , 

then no smaller frequencies need to be analyzed. 

Conversely, if, after analyzing the specified range of frequencies, any frequency in 
array co i. is equal to co r max , then a frequency larger than co r max must be analyzed in 

order to find the maximum N-Factor at that j . This can also be demonstrated by Figure 2. 
If the specified frequency range was 6000 Hz to 10000 Hz, then 10000 Hz would corre- 
spond to (0 r max . Although the N-Factor at 10000 Hz is the largest in this range, it is not 
the maximum N-Factor for this station. Larger frequencies must be analyzed. After ana- 
lyzing frequency (t> r max + Ad) r , array CD. is checked to see if co r max + Aco r appears at 

any j . If it does, then even a higher frequency must be analyzed. If it doesn’t, then N. is 

the final array of maximum N-Factors. These N-Factors form the analysis N-Factor enve- 
lope or distribution for this particular airfoil at the flow conditions. 
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2.4.2 Use of N-Factors to Estimate the Transition Location 

Since COSAL uses the e N method, the N-Factors that are calculated can be used to 
obtain an estimate of the transition location. When the N-Factors exceed a certain value, 
N tr , transition is estimated to occur at that point. 

A wide range of estimates for the value of N tr have been made, depending on whether 
the correlation was determined from a wind tunnel experiment or an in-flight experiment. 
Using wind tunnel experiments, values between nine and 1 1 have been predicted for N fr 
(refs. 23-24), although values as high as 13.5 have been estimated using in-flight experi- 
ments (refs. 10, 25). 

2.5 Airfoil Design Method 

The airfoil design method used is the CDISC (Constrained Direct Iterative Surface 
Curvature) method of Richard Campbell (ref. 26). Before designing a new airfoil, the 
method first modifies the initial target pressures to meet the desired aerodynamic and geo- 
metric constraints. Upon obtaining a final target pressure distribution, the method modi- 
fies the original airfoil to design a new airfoil that has a pressure distribution closer to the 
target pressures. 

For local Mach numbers less than 1.1, the airfoil is modified to meet the target pres- 
sures based on the relation (ref. 26) 

AC = AC p x(l + C 2 ) 05 (18) 

where C is the curvature of the airfoil, A C p is the difference between the target and analy- 
sis pressures, and t is a parameter equal to 1 for the upper surface and -1 for the lower 
surface. For local Mach numbers greater than 1.1, the relation (ref. 26) 
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( 19 ) 


-=m^x-(frr 

is initially used, where is the free-stream Mach number. As the analysis pressures 

approach the target pressures. Equation 18 is used with Equation 19 to converge the pres- 
sures more quickly. 

In addition to using Equations 18 and 19 to modify the airfoil, the angle of attack is 
adjusted based on the difference between the analysis and target pressures in the leading- 
edge region. 

Some of the parameters constrained in this airfoil design method are the lift and pitch- 
ing moment coefficients, the maximum airfoil thickness, front and rear spar thicknesses, 
the leading-edge radius, and the trailing-edge angle. 

In this application of the method, the target pressures are modified to meet the aerody- 
namic constraints prior to using CDISC. Therefore, the target pressures are modified by 
CDISC only to meet the geometric constraints, not the aerodynamic constraints. Further- 
more, the airfoil design method has been modified so that only the lower surface target 
pressures are changed to meet the geometric constraints. This has been done in an effort 
not to disturb the amount of NLF that has been achieved on the upper surface. 
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3.0 Obtaining the Target Pressure Distribution 

This chapter is devoted to the discussion of the Target Pressure Design module shown 
in the flowchart in Figure 1. This module uses the analysis N-Factor distribution from the 
stability analysis code to calculate a target N-Factor distribution that forces the boundary 
layer to transition at the desired location. Then, using the analysis and target N-Factors, as 
well as the analysis pressures, a new target pressure distribution is calculated that is closer 
to meeting the desired NLF constraints. A new airfoil can then be designed using these 
target pressures. 

The components of the Target Pressure Design module are shown in the flowchart in 
Figure 4. Each of these components will now be discussed. 

3.1 The Target N-Factor Distribution 

In order to calculate a target pressure distribution, a target N-Factor distribution that 
has the desired amount of NLF must be determined, as shown in the first module on the 
flowchart in Figure 4. 

A method has been developed for calculating a target N-Factor distribution from the 
analysis N-Factors and four control points (x , , x 2 , x cp 3 and x cp 4 ) specified in the 

streamwise direction. The first control point, x cp x , is positioned at the location where the 
analysis N-Factors first exceed an N-Factor level N cp< j . In order to calculate a realistic 
target N-Factor distribution, it is desired to retain the current analysis N-Factors as the tar- 
get N-Factors ahead of l . 

To calculate the target N-Factors aft of x cp t , the target N-Factors desired at the sec- 
ond ( N cp 2 ), third (N cp 3 ) and fourth (N cp 4 ) control points are specified. The target N- 
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Factor distribution is calculated by drawing lines through these four points and smoothing 
the curve using a polynomial fit. 

Since the shape of the target N-Factor distribution is dependent upon the speed regime 
for which the airfoil is being designed, the values for N cp x , x cp 2 , N cp 2 , x cp 3 , N cp 3 , 

x cp 4 and N cp 4 may vary from one design to the next Typical values for these parame- 
ters will now be discussed for subcritical and supercritical cases. 

3.1.1 Subcritical Cases 

Using a transition N-Factor of 10 and allowing the flow to transition at 60% chord, a 
typical target N-Factor distribution for a subcritical case is shown in Figure 5. These tar- 
get N-Factors were obtained from the analysis N-Factors shown in the same figure, the 
following control points, 

x , = 0.58 

cp, 2 

jt , = 0.60 
cp , 5 

jt , = 0.65 

cp , 4 

and the desired target N-Factors, 


at these control points. Using N y 


N . = 3 

cp , 1 


N rn , = 8 

cp, 2 


N „ , = 10 
cp, J 


N n , = 15 

cp, 4 


= 3 and the analysis N-Factors shown in the figure, 


the first control point, x cp x , is located at approximately 18% chord. 
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There are several reasons for choosing a target N-Factor distribution similar to the one 
shown in this figure. The region of the target N-Factor distribution ahead of x cp 2 forms a 

buffer zone above which the target N-Factors are not allowed to grow so that the boundary 
layer will remain laminar prior to the desired transition point at slightly off-design condi- 
tions. Beyond x cp 2 , the target N-Factors are allowed to grow rapidly to force transition. 

In this case, x cp 3 denotes the location of desired transition, which is at 60% chord for a 

transition N-Factor of 10. It is not necessary, however, for any control point to be located 
exactly at the desired transition point. 

The steep N-Factor gradient beyond x cp 2 takes into account the idea that the transi- 
tion N-Factor may not be exactly 10. In reality, the transition N-Factor could be as low as 
eight or as high as 15. Even if the transition N-Factor were eight, the airfoil in this case 
would still have NLF to 55% chord. If the transition N-Factor were 15, then the airfoil 
would have NLF to 65% chord. This indicates that the flow could actually undergo transi- 
tion anywhere between 55% and 65% chord, an uncertainty in the transition location of 
10% chord. 

Suppose that a target N-Factor distribution similar to the one in Figure 6 were used for 
a subcritical airfoil design. If the transition N-Factor is bounded between eight and 15, 
then the flow could undergo transition anywhere between 35% and 70% chord. As a 
result, this is a less desirable target N-Factor distribution since the uncertainty in the tran- 
sition location of the distribution is greater. 

There is, however, a limit on how steep the target N-Factor gradient beyond x 2 can 
be. If the N-Factor gradient is too great, then the adverse pressure gradient required to 
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obtain the target N-Factors in that region will cause the laminar boundary layer to separate 
in that region. Experience has shown that in order to avoid laminar separation the target 
N-Factors should not be allowed to grow more than 10 N-Factors for every 10% chord. 
3.1.2 Supercritical Cases 

In supercritical designs, it is much more difficult to find a realistic target N-Factor dis- 
tribution. The target N-Factor distribution shown in Figure 5 is not realistic for supercriti- 
cal airfoil designs. Since the large N-Factor gradient aft of x cp 2 would cause an adverse 
pressure gradient in the target pressure distribution, it is likely that a shock may form 
ahead of x 4 . In the case of a supercritical airfoil design, this is unwanted since x cp 4 


denotes the desired location of the shock. 

Target N-Factor distributions similar to that seen in Figure 6 were initially used in the 
design of airfoils in this supercritical regime; however, using these target N-Factor distri- 
butions resulted in unrealistic target pressure distributions, which caused a shock to form 
ahead of x cp 4 . As a result, since the analysis pressures would never converge to the tar- 


get pressures in the design of an airfoil, the analysis N-Factors would never converge to 
the target N-Factors. 

After running many cases, a target N-Factor distribution similar to the one shown in 
Figure 7 was found to be a realistic distribution for supercritical cases. This distribution 
was obtained using the following control points 

x , = 0.10 

x , = 0.35 

cp, 3 


x 4 = 0.55 

cp, 4 


19 



The location of x .is at 0% chord since jV , = -1 and the analysis N-Factor at 0% 

CPj 1 * 

chord is 0. In this figure, the target N-Factors are allowed to grow rapidly from 10% to 
35% chord, where the flow can tolerate the mild adverse pressure gradient that is required 
to obtain this N-Factor distribution. Then, beyond 35% chord the N-Factors are allowed 
to grow only by a small amount so that a shock will not form ahead of x cp 4 . 

A result of using a target N-Factor distribution similar to the one shown in Figure 7 is 
that the uncertainty in the transition location can be large. Since a steep N-Factor gradient 
is not realistic in supercritical designs, an N-Factor distribution that accounts for the 
uncertainty in the transition N-Factor is not possible since it is difficult to force the N-Fac- 
tors to grow much higher than eight or nine. As a result, the design philosophy for super- 
critical cases is to assume a transition N-Factor near eight. In doing this, the boundary 
layer will remain laminar until the flow encounters a shock at x cp 4 . An unfortunate result 

of this idea is that if the transition N-Factor is 13 in reality, then laminar separation will 
occur at the shock. 

Note how there are no large N-Factor gradients in the distribution beyond 35% chord, 
which would indicate that there should not be any large changes in the pressure gradient 



on the resulting airfoil. Thus, the airfoil that results should have a flat, roof-top pressure 
distribution that is typical of many supercritical airfoils. 

3.2 Extrapolation of Analysis N-Factors 

The laminar boundary layer solver is valid only for attached, laminar boundary layers, 
and is terminated once the laminar boundary layer separates. As a result, the boundary 
layer solver only supplies data to the stability analysis code as long as the flow is attached. 
This in turn means that the stability analysis code can only calculate N-Factors until the 
boundary layer separates. 

Suppose that the NLF constraint requires that the flow be attached to 60% chord, but 
the current airfoil’s boundary layer separates at 30% chord. This means that the stability 
analysis code can only calculate N-Factors to 30% chord, although a target N-Factor dis- 
tribution will be specified to 60% chord. The N-Factor design method requires analysis 
and target N-Factors at every station ahead of the fourth control point, which is located at 
60% chord in this case. This indicates that the analysis N-Factors must be extended 
through the separated flow so that the N-Factor design method discussed in Section 3.3.1 
can be used. As a result, the next module in the flowchart shown in Figure 4 is used for 
this purpose. 

The method that is used to perform the extrapolation depends on whether the flow is 
subcritical or supercritical. For subcritical cases, this is accomplished by linear extrapola- 
tion based on the analysis pressure distribution. With station j = k being the current 
location of laminar separation and j = l being the location of the fourth control point, the 
N-Factors between k + 1 and l can be calculated as 
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( 20 ) 


Ns _ < — N j _ 2 

Nj - Nj_ \ + Q _ Q ( Cp a j u - Cp a j_l u ) 

p, a,j - 1 , u p,a,j-2,u 


where Nj denotes the current analysis N-Factor distribution and C p a j u denotes the cur- 
rent analysis pressure distribution. From experience, for stability of the N-Factor design 
method the N-Factors calculated by Equation 20 are restricted as follows: 


+ 150(Xj-Xj_ 1 ) 

( 21 ) 

Nj> 0 

( 22 ) 


For supercritical cases, a different approach is taken. In supercritical flows, it is 
unlikely that the laminar boundary layer will remain attached through the shock. More- 
over, a shock wave is a discontinuity in the pressures and would cause difficulties for the 
N-Factor design method that is used to design a target pressure distribution. Therefore, 
calculating N-Factors aft of a shock wave, even if the flow remained attached through the 
shock, is unnecessary. As a result, in the case of supercritical designs, the laminar bound- 
ary layer solver is terminated at the station upstream of the shock wave. 

In the case where a shock wave exists ahead of the fourth control point, it is necessary 
to move the shock aft to the desired location, which is at the fourth control point, so that 
the N-Factor design method can be used. This is done by imposing a flat roof-top pressure 
distribution between the current shock location and the fourth control point. If the current 
shock is at station j = k and the new shock is to be located at j = l (the station repre- 
senting the fourth control point), then the current shock’s upstream pressure coefficient is 
maintained aft to station / . That is, the analysis pressures between k + 1 and / are rewrit- 
ten as 


r* + 1 - r' 

p,a,j,u p, a, k, u 


(23) 


22 



where u 3X6 new analysis pressures and C? p a j u are the original analysis pres- 

sures. This is demonstrated in Figure 8. 

Also in this case, analysis N-Factors only exist ahead of station k since the laminar 
boundary layer solver is terminated at k due to the shock wave. This indicates a need to 
extrapolate the analysis N-Factors aft to station / so that the N-Factor design method can 
be used to calculate a target pressure distribution. Since the pressure coefficients between 
k + 1 and / are the same, the extrapolation scheme described in Equation 20 cannot be 
used to calculate a new analysis N-Factor distribution. Therefore, it is assumed that the N- 
Factors grow only with x , and hence can be written as 


N. . -N. 2 

N f = N. l + } — (jc - - Xj _ x ) (24) 

J 1 1 X j-l- x j-2 J J 
between stations k + l and l . 

It is possible, however, that the laminar boundary layer will separate before the shock. 
In this case, prior to using Equation 24, Equations 20-22 must be used between stations 
m + 1 and k , where m is the location of laminar separation and k is the location of the 


shock. 


3.3 The Upper Surface Target Pressures 

Now that the analysis and target N-Factor distributions are known, a target pressure 

distribution can be calculated so that the CDISC airfoil design method can be used to 


design a new airfoil. The first step in this process is to calculate the upper surface target 
pressure distribution. In order to do this, the upper surface pressures are divided into two 
important regions, as shown in Figure 9. In the first region, a new N-Factor design 
method is used to calculate the target pressures based on the difference between the target 
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and analysis N-Factors. The pressures in the second region are calculated from the recov- 
ery pressures of the first airfoil that was analyzed. The specifics of how the pressures are 
determined in these two regions will now be discussed. 

33.1 The N-Factor Design Method 

The next module on the flowchart of Figure 4 is used to calculate the upper surface tar- 
get pressures in the region that extends from the stagnation point of the airfoil to the fourth 
control point. This is the region where the current pressure coefficients are modified to 
move the current N-Factor distribution towards the target N-Factor distribution that was 
prescribed earlier. In order to do this, a new method had to be developed. 

In order to achieve the desired amount of natural laminar flow, a method was devel- 
oped that possessed the following properties: 

1. If the current N-Factor at a given airfoil station was larger/smaller than the 
target N-Factor at that station, then the pressure gradient would have to 
become more favorable/adverse at that station in order to decrease/increase 
the current N-Factor. 

2. The N-Factor at a given point on the airfoil would be changed by modify- 
ing only the pressure coefficient at that point. (This local change in pres- 
sure coefficient would have an effect on the N-Factors downstream of the 
current airfoil station, but not the N-Factors upstream since the boundary 
layer equations are parabolic.) 

3. The design method must produce a smooth and continuous pressure distri- 
bution between the stagnation point and the fourth control point. 
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With these ideas in mind, this method, called the N-Factor design method, was devel- 
oped such that the change in pressure coefficient required at a given airfoil station j would 
be governed by the linear relation 

AC . 

/'*' “ = A AN- (25) 

Ax -j J 

where 


AC • = C t - C 

p,j,u p, T,j, u p,a,j,u 


Axj = \x rXj _ x \ 


ANj = N T j-Nj 


(26) 

(27) 

(28) 


In these equations, A is a relaxation factor (typically 0.50), C p T . u are the new upper 
surface target pressures, C are the current upper surface analysis pressures. Ax- is 
the grid spacing, N T . are the target N-Factors, and Nj arc the current analysis N-Factors. 

Note that Equation 25 satisfies condition 1 above in that for a positive A Nj, A C p ■ u is 
also positive, which induces a more adverse pressure gradient at station j . Conversely, for 


a negative AN } , A C p ■ u is also negative, which produces a more favorable pressure gra- 


dient at station j . 

Now in order to satisfy condition 2 above, the N-Factors downstream of station j must 
be corrected for the change in pressure that has been imposed. This is accomplished by 
assuming that each N-Factor downstream of station j is increased or decreased by as 
much as the N-Factor at station j was changed. Effectively, then, ANj is added to each N- 

Factor downstream of j . Condition 3 above requires that the pressure distribution be both 
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smooth and continuous. In order to maintain a continuous pressure distribution, the 
change in pressure applied to station j is also applied to each station downstream of j . In 
doing this, the pressure distribution remains inherently smooth. 

Condition 2 indicates that each N-Factor downstream of j must be corrected for the 
change in pressure coefficient at j . This perspective can be reversed. If the method is 
designing at station j , all the modifications in pressures and N-Factors upstream must be 
taken into account at station j. By solving Equation 25 for C p T j - u and correcting for 
upstream effects as mentioned above, the new relation 


is obtained. With j = 1 corresponding to the stagnation point, this relation is valid from 


j = 2 to j = / , the location of the fourth control point. The boundary conditions at the 
stagnation point are 

Cp,T, l,u = ^p,a,l,u = Cp' Sfag (30) 

N t ;1 = N t = 0 (31) 

Using COSAL as the stability analysis code, Equation 31 is automatically satisfied as long 
as the target N-Factor at the stagnation point is set equal to zero. 

Figure 10 illustrates how the N-Factor design method is used to calculate the new tar- 
get pressures. The left-most plot in Figure 10 shows sample analysis and target N-Factor 
distributions. The change in N-Factor, AN , that is required at a grid point is calculated as 
the difference between the target and analysis N-Factors at that point. The A C p that is 

necessary to make this change in N-Factor is then calculated using Equation 25. Then, as 
shown in the center plot in Figure 10, the pressure coefficient at the current grid point and 
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each grid point downstream is changed by A C p . The analysis N-Factors at the current 

grid point and each one downstream are changed by AN to obtain a modified analysis N- 
Factor distribution. This is illustrated in the right-most plot of Figure 10. These plots 
demonstrate how the N-Factor design method is used at one grid point. The same 
approach is used at each of the grid points over which a target N-Factor distribution is 
specified. 

As mentioned above, the parameter A used in Equation 25 is a relaxation factor in the 
design process. Using COSAL as the stability analysis code, a typical value for A is 0.50. 
Although it is not recommended for stability purposes that A be increased above 0.70, A 
may be decreased if more stability is desired. In addition, the value of A does not seem to 
be affected by flow parameters like Mach number or Reynolds number. 

Perhaps some physical significance can be found in the relationship of Equation 25. 
As the Tollmien-Schlichting waves propagate downstream, their velocity is 


v — c + U , 

W W ’ 


(32) 


where c w is the phase velocity of the waves and U is the local velocity of the flow. If the 

flow is being accelerated (i.e., the velocity U is increasing), then the speed of the wave, 
v w , is also increasing. As v w increases, the T-S waves spread out and their amplitudes 
decrease, which promotes boundary layer stability. On the contrary, if the flow is being 
decelerated (i.e., the velocity U is decreasing), then the speed of the wave is also decreas- 
ing. As v w decreases, the T-S waves bunch-up and their amplitudes increase, which 

increases the instability of the boundary layer. This indicates that there is a negative pro- 
portional relationship between a change in the amplitude of the T-S waves and a change in 
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the local flow velocity. Since the N-Factor is merely the logarithmic growth of the ampli- 
tude of the T-S waves, then there is also a negative proportional relationship between a 
change in N-Factor and a change in the local flow velocity. 

There is also a negative inverse relationship between a change in the local pressure 
gradient and a change in the local flow velocity. As a result, there is a direct relationship 
between a change in the local pressure gradient and a change in N-Factor, as proposed in 
Equation 25. 

3.3.2 The Pressure Recovery Region 

The next module on the flowchart of Figure 4 is used to calculate the upper surface tar- 
get pressures in the recovery region. The pressure recovery region is composed of the 
pressures between station j = / and the trailing-edge of the airfoil (J = n ), as shown in 
Figure 9. These pressures are formed by modifying the recovery pressures of the first air- 
foil that was analyzed. C . 0 denotes the upper surface pressure coefficients of the first 

airfoil that was analyzed, with j denoting the airfoil station. From this pressure recovery, 
two intermediate pressure distributions are formed and used to determine the final target 
pressures in the recovery region. Figure 11 shows this process. 

The first intermediate recovery pressure distribution is determined by linearly scaling 
(see Appendix A, page 89) C p ; 0 , as shown in Figure 1 1 (a). The target pressure coeffi- 
cient at j — l was already determined by the N-Factor design method described in the 
previous section. Therefore, C p . 0 must be linearly scaled so that its pressure coefficient 

at station / is the same as the pressure coefficient at station / given by the N-Factor design 
method. As a result, the first intermediate recovery distribution is given by 
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(33) 


"p.y. i 


C , , ~C „ 

_ p, 7, /, u p, n, 0 

r _ r 

u p, /, 0 °p, n, 0 


( Cp f j, 0 ~ Cp, n. o) + Cp, n, 0 


Equation 33 is valid from j = l to j = n . 

The second intermediate recovery distribution is obtained by adding a linear loading 
distribution onto C p j 0 , as shown in Figure 11 (b). This loading distribution is added in 

order to match the design pitching moment coefficient constraint. The process by which 
this is done will be described in Section 3.4.2. The second intermediate recovery distribu- 
tion is then given as 


r - r + - 
pJ< 2 pJ'O 0.30 


1-Jt. 

J C 


p, 0.7c 


(34) 


where C p 0 7c denotes the magnitude of the loading distribution at 70% chord. This equa- 


tion is also valid from j = l to j = n. 

Now, the final target pressures of the recovery region are obtained by taking a 
weighted average (see Appendix B, page 91) of the first and second intermediate recover- 
ies, as shown in Figure 1 1 (c). In doing so, the final target pressures are given by 


x n -x- 

C T - =— l C ■,+ J - -C 

p.T.J.U X -X. P’J - 1 X -X. P-J 


X - - X, 


X n~ X l 


X n~ X l 


(35) 


Being valid from j = / to j = n , this expression allows C p T ■ u to retain the value of 


C p j j at j = l and the value of C p j 2 at j = n. 

The procedure described in this section is only valid for subcritical flows. For cases 
where there is a shock wave on the upper surface (at j = / by design), the same process is 
used to determine the pressure recovery region except that everywhere an “/” appears in 
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the equations, “/ + 1 ” should be put in its place. As a result, then, the equations are only 
valid from ; = / + 1 to j = n for supercritical cases. 

In order to use these equations from j = / + 1 to j = n , the target pressure coefficient 
on the downstream side of the shock (J = l+l) must be calculated using the free-stream 
Mach number . The relation (ref. 27) 


2 f 1 +0.2mJY 5 

C = — - — — - 1 

P lAM 2 ^ l +0.2M 2 j 



2 ffl +0.2M 2 Y*' 5 

C = — - — — - 1 

P ’ 2 1.4M 2 U 1 +0.2MI J J 

oo Z 


(39) 


However, because of curvature of the airfoil, the pressure coefficient calculated by 


Equation 39 is not totally correct. To account for curvature effects, the target pressure 
coefficient on the downstream side of the shock is given as 
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( 40 ) 


Cp,T,l+l,u - 2^p* + Cp,2) 

where C p * is the sonic pressure coefficient. The process previously discussed can now be 

used to determine the remainder of the target pressures in the recovery region. 

3.4 Meeting the Aerodynamic Constraints 

Once this preliminary target pressure distribution is determined, the next module on 

the flowchart of Figure 4 is used to modify the upper surface target pressures to meet the 
aerodynamic constraints, which are the lift and pitching moment coefficients. A process 
was developed by which the upper surface target pressures could be modified to meet 
these constraints, while changing the N-Factor distribution as little as possible so that the 
NLF would not be disturbed. 

To match these constraints, the upper surface target pressure distribution is divided 
into three segments. Figure 12 shows a typical upper surface target pressure and target N- 
Factor distribution, with the pressure distribution being divided into the following three 
regions: 

1. The leading-edge region. This region extends from the stagnation point 0 = 1 ) to the 
last point where the target N-Factors are zero (j = m ). 

2. The center region. This region extends from airfoil station m to the fourth control 

point 0‘ = /)• 

3. The pressure recovery region. This region extends from station / to the trailing-edge 
0 = n). 
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The methods by which the target pressures are modified to meet these aerodynamic 
constraints are discussed in the next two sections. A flowchart showing the procedures of 
the sections appears as Figure 13. 

3.4.1 Lift Coefficient 

Since the extent of NLF is dependent upon the pressure gradient, a method for modify- 
ing the upper surface target pressures to match the lift coefficient was developed that 
would maintain the current pressure gradient through the region where the N-Factors are 
increasing most rapidly. The center region is the region where the N-Factors are most 
important because they are growing fastest. As a result, it is desired to shift each pressure 
coefficient in this region by the same amount in order to preserve the current pressure gra- 
dient. But, if these pressures are shifted, then the pressures in the leading-edge and pres- 
sure recovery regions must also be manipulated so that the pressure distribution remains 
smooth. In the leading-edge region, the pressures are linearly scaled so that they are con- 
tinuous at j = m , the beginning of the center region. 

In order to determine how much the current target pressures C p T - u are to be modi- 


fied, it is first assumed that the center region is going to be shifted by A C p c . This would 


cause a change in lift coefficient of 

Ac. = AC , 

/, C p,C i 


S l~ S m 


( 41 ) 


V “n J 

where s. represents the chordwise distance from the stagnation point 0 = 1) defined as 


S- 


7 = 

* = 2 


( 42 ) 
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If the pressures in the center region are shifted by A C p c , then the pressures in the lead- 


ing-edge and pressure recovery regions must also be changed, as will be discussed 
momentarily. The change in lift coefficient in the leading-edge and pressure recovery 
regions may be approximated by a linear loading variation along the corresponding arc 
length. Therefore, 


( . \ 


Ac, . = -AC 

l, le 2 p, c 


\ S nJ 


(43) 


Ac, nr = rAC 
Ipr 2 P' c 


• c _ c ' 

S n S l 


V " J 


(44) 


As a result, the total change in Ac l for these changes in pressure would be the sum of 


Equations 41, 43, and 44, which results in 


Ac / = AC P , C | 


Fm . S l~ S m . l S n~ S l 


2s s 

V n 


+ - 


2 S 


(45) 


n J 


Simplifying and solving Equation 45 for AC^ c , the amount that the target pressures in 
the center pressure region must be shifted by is 


AC 


2r Ac, 

n l 


p ’ c s, s + s 

i fu n 


(46) 


For stability reasons, this equation can be rewritten to include a relaxation factor. In addi- 
tion, the change in lift coefficient that is needed is actually the difference between the 
design lift coefficient and the lift coefficient of the current target pressures. So, with a 
relaxation factor X (typically 50%), Equation 46 becomes 


AC 

p,c 


C /) 

S l~ S m + S n 


(47) 
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Thus, after calculating A C p c , the new target pressures over the upper surface from 


j = m to j = l become 


C* 


+ 1 

p, T,j, u 


= c 


p. T.j, 


+ AC 

U p, C 


From j = 1 to y = m - 1 , the upper surface target pressures become 


(48) 


C l + 1 _ (f f . ■ \ ■ 

^J+l _ ^p,T,m,u '■'p, T, \,u[ ^ 1 , r* 

l p,tj,u ~ * _i \^p.tj,u Vu ^ S . r . i .* 


(49) 


C‘ - c 

'■'p, T, m, a '•'p, T, 1 , u 


The region between j = l + 1 and j = n is the pressure recovery region. Therefore, the 
method originally used to obtain the pressures in the recovery region (described in 


Section 3.3.2) is again used to obtain 


d V . 

P,T,j,u 


from j = l to j = n . 


As shown in Figure 13, this process is then repeated until 


C l, des C l ^ C l, tol 


(50) 


where c t tol is the desired tolerance for the lift coefficient (typically 0.01). 

3.4.2 Pitching Moment Coefficient 

Once the design lift coefficient has been achieved, a modification to the pressures in 
the recovery region is made in an attempt to match the design pitching moment coeffi- 
cient. The method by which the target pressures in the recovery region are determined is 
described in Section 3.3.2. In that section, a linear loading distribution with magnitude 
C p o7c at 70% chord was included in Equation 34 for the second intermediate pressure 


recovery. This exists for the sole purpose of achieving the design pitching moment coeffi- 
cient. 

In order to meet the pitching moment constraint, lift is transferred to or from the pres- 
sure recovery region by changing the amount of lift in the loading distribution. This is 
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accomplished by modifying C p Qlc by adding A C p 07c . The change in lift coefficient of 
the target pressures resulting from adding A C p 07c to C p 07c is 

A c i = ~2 n ~ x i) o.7c (51) 


where x { represents the location of the fourth control point. This change in lift coefficient 
causes a change in pitching moment coefficient 

Ac m = -Ac z (x z - 0.25 + (52) 

Using Equation 51, Equation 52 can be simplified to 

Ac m = |( 1 -^)(^ / + ^AC p0 .7c < 53 > 

Simplifying and solving for A C p Qlc , 


AC 


24Ac_ 


P’ Qlc ( 1 -*,) ( 8 jc z + 1 ) 


(54) 


Since the change in pitching moment required is actually the difference between c m des 


and c m , A c m can be replaced with c m dgs - c m . In addition. Equation 54 can be written to 


include a relaxation factor q (typically 50%). As a result, the final expression for 

0.7c 


AC 


P’ 0Jc ( 1 — X[) (8x z + 1) 


(55) 


Therefore, the new value of C p 0 lc is 


C p,0.7c “ t'p, 0.7c + ^p. 0.7c 


(56) 
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where Cf 0 7c is the original value. Using this new value of C p 01c , the target pressures 

in the recovery region can be recalculated using the method of Section 3.3.2. 

After changing the pressures in the recovery region to move c m towards c m des , the 

lift coefficient of the target pressure distribution was changed. As a result, the method of 
Section 3.4. 1 is used to modify the target pressures to once again achieve the design lift 
coefficient. 

This process is repeated until 


^ l, des 


~ C. < C 


l, lol 


(57) 


and 


c . — c ^ c . , 

m, des m m, tol 


(58) 


A typical value for c m {ol is 0.01. 

3.5 Adjusting the Leading-edge Radius 

Since the airfoil design method uses a target pressure distribution to design a new air- 
foil, the leading-edge radius of the redesigned airfoil depends on the shape of both the 
upper and lower surface target pressure distributions. As a result, many different target 
pressure distributions can be used to meet the leading-edge radius constraint. However, it 
is desirable to use a target pressure distribution that has a reasonable amount of NLF on 
the lower surface, since the N-Factor design method described in Section 3.3. 1 is not used 
to design the lower surface target pressures to obtain NLF. 

In order to do this, the pressure distribution shown in Figure 14 is used to increase the 
leading-edge radius of the airfoil through modifying the leading-edge target pressures of 
the upper surface. These pressures were the leading-edge pressures of an airfoil that was 
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redesigned starting from a NACA 64 r 212 airfoil. The redesigned airfoil had a large 
extent of NLF on the lower surface and these upper surface pressures seemed to be a con- 
tributor to this. 

To use these leading-edge pressures, this distribution is linearly scaled so that its pres- 
sure coefficient at 30% chord is the same as that of the current target pressures. Then, 
over the first 30% of the chord, a weighted averaging technique is used to calculate the 
new target pressures using this distribution and the current target pressures. The pressure 
coefficient at the leading-edge of the distribution shown in Figure 14 is retained as the new 
target pressure, while the current target pressure coefficient at 30% chord remains 
unchanged. 

After modifying the target pressures through this averaging, the method of Section 3.4 
must be used to modify the upper surface target pressures to meet the aerodynamic con- 
straints. This is shown in the next module of the flowchart of Figure 4. 

3.6 The Lower Surface Target Pressures 

The next module on the flowchart of Figure 4 is used to calculate the lower surface tar- 
get pressure distribution. Although the N-Factor design method is not applied to design 
the lower surface since the lower surface target pressures are modified to meet the geomet- 
ric constraints, it is desired to obtain as much NLF on the lower surface as possible. In 
order to do this, the upper surface target pressure distribution is linearly scaled and a 
weighted average is taken in order to obtain the initial lower surface target pressures. The 
following process is only applied on the first iteration of the method since it has been seen 
occasionally to cause the target pressure distribution to change too much from one itera- 
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tion to the next when many constraints have been imposed. Figure 15 demonstrates the 
process that is about to be discussed. 

The intermediate lower surface target pressures are obtained by linearly scaling the 
upper surface target pressures, as demonstrated in Figure 15 (a). If C p T j u represents the 

upper surface target pressures, then the intermediate lower surface target pressures are 
given as 


,, _ Cp'T, 1 ■ K ~ Cp, T, k, u k C l, des 

'p,j, 1 C T , - C T t 

p, 7, l,u p, T, k , u 


( C p , r,j, u~C p T i „) + C p T j u (59) 


where k denotes the airfoil station closest to 50% chord, and A C p k is used to cause the 


final target pressure distribution to have a lift coefficient of c ; des . Equation 59 is valid 


from y=l (the stagnation point) through j = n (the trailing-edge). 

For the region on the lower surface forward of the fourth control point on the upper 
surface, the pressures represented by Equation 59 are used for the target pressures in the 
recovery region. That is, from j = 1 to j - l (the location of the fourth control point), 
the lower surface target pressures are given by 


C Tl = C , 

p,T,j,l p,),\ 


(60) 


For the region aft of station / , the final lower surface target pressures are found by taking 
a weighted average of the intermediate pressures and the analysis pressures of the recov- 
ery region of the initial airfoil that was analyzed, which are denoted by C p ■ 0 . This is 


demonstrated in Figure 15 (b). As a result, from j = / + 1 to j = n t , the lower surface 
target pressures become 
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(61) 


-C_ ; 


P* T,j, l s _c p,j,0 s p,j, 1 

ft. L n, i 


Equation 61 allows C p T • t to retain the value of C p ■ x at j = / and the value of C p ^ 0 

at the trailing-edge. The final target pressure distribution is shown in Figure 15 (c). 

It should be mentioned here that these lower surface target pressures will be modified 
to meet the geometric constraints. Therefore, the amount of NLF that is obtained on the 
lower surface is directly a function of the geometric constraints imposed. Moreover, since 
the geometric and aerodynamic constraints constantly react to the changes that the other 
makes to the pressures, the amount of NLF obtained on the lower surface is also depen- 


dent upon the aerodynamic and upper surface NLF constraints. 
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4.0 Modifying Target Pressures to Enforce Constraints 

The flowchart of Figure 1 shows a module above the airfoil design method that is 

labeled Modify Target Pressures to Enforce Constraints. After a target pressure distribu- 
tion is calculated as described in Chapter 3, these taiget pressures are modified while a 
new airfoil is being designed by the airfoil design method. These changes are needed in 
order to enforce the desired aerodynamic and geometric constraints. 

A flowchart of the Modify Target Pressures to Enforce Constraints module is shown in 
Figure 16. Each of the components in this module will now be discussed. 

4.1 Leading-Edge Pressures 

The CDISC airfoil design method directly modifies the airfoil to meet the leading- 
edge radius constraint, without modifying the target pressure distribution as it does to 
meet other constraints. As a result, the target pressures are modified within the first mod- 
ule of Figure 16 to account for the change that the airfoil design method made to the lead- 
ing-edge of the airfoil to meet the leading-edge radius constraint. 

This module is not used until 10 iterations through the airfoil design method have been 
completed. This allows the analysis pressures to approach the desired target pressures in 
the leading-edge region. But, without modifying the target pressures, the analysis pres- 
sures would never exactly match the target pressures in this region, while still maintaining 
the leading-edge radius constraint. 

After completing the first 10 iterations, only the lower surface target pressures are 
modified. This is done by taking a weighted average ahead of 30% chord of the current 
lower surface analysis pressures and the current lower surface target pressures. As a 
result, the new lower surface target pressures become 
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where k represents the station nearest to 30% chord. 

During the second 10 iterations through the airfoil design method, only the lower sur- 
face target pressures are modified. This allows the upper surface analysis pressures to 
approach the upper surface target pressures, which have been designed to meet the NLF 
constraints. After these first 20 iterations have been completed, the upper surface analysis 
pressures should closely resemble the upper surface target pressures. The new upper sur- 
face target pressures are then calculated using a weighted average between the current tar- 
get and analysis pressures. Therefore, the new upper surface target pressures can be 
calculated as 


c' V • = i - - 1 c 1 

p,T,j,u s p, a,j, u 

\ k) 




where k once again denotes the station nearest to 30% chord. In doing this, the upper sur- 
face target pressures will usually not change much since the airfoil design method has 
completed 20 iterations. 

4.2 Trailing-Edge Pressures 

To meet the trailing-edge angle constraint, the airfoil is modified directly, without 
changing the target pressures. Therefore, in the next module shown in Figure 16, the trail- 
ing-edge target pressures also need to be modified to account for the changes that have 
been made to the airfoil to meet the trailing-edge angle constraint. On the lower surface, 
the new target pressures are calculated from the current target and analysis pressures using 
a weighted average aft of 60% chord. On the upper surface, the new target pressures are 
calculated from the current target and analysis pressures using a weighted averaged tech- 
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nique aft of the fourth control point, which is at station /. In the form of equations, the 
new upper and lower surface target pressures become 


+ 1 

'P.TJJ 


S j ~ S _k (4 

s n~ s k P ' 
n i k 


S. -5 




a ’j ' 1 s n S-P-T.JJ 


+ 1 


Sj V 


s n~ s j r j 

' P’TJ.u ~ s _ _ s ^p,a,j, u + S "-s t P.T.j.u 


(64) 


(65) 


'/ 


where k is the station nearest 60% chord on the lower surface, and / is the location of the 
fourth control point on the upper surface. These modifications are made at the same fre- 
quency that the leading-edge target pressures are modified to account for the leading-edge 
radius constraint. 

The next module on the flowchart of Figure 16 is used to modify the target pressures to 
meet the aerodynamic constraints. This process was discussed in Section 3.4. 


4.3 Releasing Constraints 

After modifying the target pressures using the methods of the previous two sections, 
the next module on the flowchart is used to release one or more of the constraints in the 
event that the design is over-constrained. In order to do this, a constraint priority list was 
established. It is as follows: 

1. Upper surface NLF 

2. Section lift coefficient 

3. Leading-edge radius 

4. Spar thicknesses 

5. Maximum airfoil thickness 

6. Trailing-edge angle 
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7. Pitching moment coefficient 

The pitching moment coefficient is the least important constraint and would be the first to 
be released. 

Throughout the design process, the target pressures are constantly being modified to 
meet each of the design constraints. If the design is not over-constrained, the target pres- 
sures will move toward a distribution that satisfies all of the constraints. It is possible, 
however, that the problem is over-constrained and the target pressures do not move toward 
any single distribution. It is in this case that one or more of the constraints must be 
released to allow the target pressures to approach a distribution that meets the more impor- 
tant design constraints. 

In order to determine if any of the constraints need to be released, the average amount 
that the target pressures change by from one iteration to the next is calculated. This 
change is given by the relation 


AC = Y Ic* r . - C‘ 1 ■ + y T 1 1 ~ C, ' r i J / (n + n, - 1 ) (66) 

p, av I p,T,j,u p, T,j, u ZZ p, T,j, l P,T,J,1 l 


where i denotes the current iteration. If, after 30 iterations through the CDISC airfoil 
design method, 

AC^O.IO, (67) 

then a constraint must be released. This criterion was established after running many 


design cases. It was found that when the design was over-constrained, A C p av was con- 
sistently near 0.30 or 0.40 since the target pressures were constantly changing. In the 
cases where the design was not over-constrained and the target pressures approached a 

single distribution, A C p av was usually 0.001 or smaller. 
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There arc two reasons why a constraint may be released. First, it just may not be pos- 
sible to maintain all of the constraints within the specified tolerances, which are typically 


C ltol : 

= 0.01 

(68) 

tol 

= 0.01 

(69) 

l lol ~ 

0.001 

(70) 

le, tol “ 

: 0.0001 

(71) 

tol 

= 0.50 

(72) 


The cases that are most sensitive to this are the incompressible subcritical cases and the 
supercritical cases. In these cases, even using different target N-Factor distributions can 
cause various constraints to be released. This is even more likely when the specified target 
N-Factor distribution is not realistic. 

Next, constraints may be released merely because the design constraints imposed on 
the airfoil may not be realistic. If one wishes to design a 15% thick airfoil for general avi- 
ation applications, but also wants the airfoil to be 9% thick at 20% chord and 7% thick at 
60% chord, it is very possible that the thickness constraint is not realistic for the specified 
front and rear spar constraints. In this case, the thickness constraint would be released 
according to the constraint priority list specified above. 
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5.0 Results of the NLF Airfoil Design Method 

The NLF airfoil design method described in the previous chapters has been used to 
design many airfoils. Airfoils have been designed for Mach numbers between 0.10 and 
0.80, Reynolds numbers between 1 and 20 million, and maximum thicknesses between 
10% and 18% chord. Although the majority of these airfoils have been designed for sub- 
critical cases, a few have been designed for supercritical cases where a shock wave exists 
on the upper surface of the airfoil. 

To demonstrate the method, the results of airfoils designed for a glider, commuter and 
subsonic transport aircraft are presented in the following three sections. 

5.1 Airfoil for a Glider 

For the first design case, the NACA 64 j -2 12 airfoil was redesigned at the following 
flow conditions: 

AC = 0.10 
Re = 3 million 

These flow conditions are typical of a glider. For these conditions, the pressure distribu- 
tion of the NACA 64 r 212 airfoil is shown in Figure 17, while the upper and lower N-Fac- 
tor distributions are shown in Figure 18. With these flow conditions and a transition N- 
Factor of 13.5, the laminar boundary layer for both surfaces of this airfoil separated 
instead of undergoing transition. Since laminar separation is not desired in the design of 
an airfoil, in the process of redesigning this airfoil a target N-Factor distribution was spec- 
ified that would force boundary layer transition rather than allow the laminar boundary 
layer to separate. 
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It is desired that the redesigned airfoil have the following aerodynamic and geometric 
characteristics: 

x, ru = 0.65 
C[ = 0.30 
c m = -0.060 
t max - 0.150 
t = 0.120 at x = 0.20 
t = 0.090 at x = 0.70 
r k = 0.0140 

Using the NLF airfoil design method, the upper surface analysis N-Factor distribution 
shown in Figure 18 and the target N-Factor distribution shown in Figure 19 were used to 
calculate the pressure distribution shown in Figure 20. Then, using the airfoil design 
method, the airfoil shown in Figure 20 was designed to meet the target pressures calcu- 
lated by the NLF airfoil design method. 

The redesigned airfoil is compared with the NACA 64 r 212 airfoil in Figure 21. In 
addition, some of the characteristics of the NACA 64 r 212 and the redesigned airfoil are 
compared in Table 1. In the design of this airfoil, the pitching moment and maximum air- 
foil thickness constraints were released by the process described in Section 4.3. Neverthe- 
less, the design pitching moment coefficient was coincidentally achieved. The fact that 
the maximum airfoil thickness constraint was released implies that it was probably not a 
realistic constraint given the desired front and rear spar thickness constraints. 

Only five iterations of the method were required to design the new airfoil, which took 
nearly four hours on a Silicon Graphics Indigo2 workstation with an R4000 processor. 
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Table 1. A comparison of the design constraints and the 
characteristics of the NACA 64] -2 12 and the redesigned airfoil 
at AC = 0.10, Re = 3 million, and c/ = 0.30 



NACA 64]-212 

REDESIGN 

CONSTRAINTS 

a 

1.48 

0.67 


c m 

-0.034 

-0.062 

-0.060 

Cd 

0.0044 

0.0032 ; 


x tr,u 

0.48 (lam. sep.) 

0.66 

0.65 

x tr.l 

0.53 (lam. sep.) 

0.51 (lam. sep.) 


^max 

0.120 

0.142 

0.150 

t at x = 0.20 

0.104 

0.119 

0.120 

t at x = 0.70 

0.067 

0.090 

0.090 

r le 

0.0108 

0.0140 

0.0140 


5.2 Airfoil for a Commuter Aircraft 

As the next example, the NLF airfoil design method was used to redesign the NACA 
1412 airfoil at a subcritical speed, with the flow conditions being 

AC = 0.60 
Re = 20 million 

These flow conditions are representative of a commuter aircraft. 

The new airfoil was to have the following design characteristics: 

x tnu = 0.60 
c, = 0.40 
c m = -0.080 
W® 0.120 
t = 0.095 at x = 0.20 
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t = 0.070 at x = 0.70 


= 0.0100 

For these flow conditions, the pressure distribution of the NACA 1412 airfoil at the design 
lift coefficient is shown in Figure 22. In addition, the upper and lower surface N-Factor 
distributions for this airfoil at these conditions are shown in Figure 23. In this figure, an 
N-Factor of 10 is used to determine where transition from laminar to turbulent flow occurs 
in the boundary layer. 

After calculating the target N-Factor distribution shown in Figure 24, the NLF airfoil 
design method calculated the target pressure distribution shown in Figure 25. Using this 
pressure distribution, the airfoil shown in Figure 25 was designed by the CDISC airfoil 
design method. Figure 26 shows a comparison of the NACA 1412 airfoil and the rede- 
signed airfoil. 


Table 2. A comparison of the design constraints and the 
characteristics of the NACA 1412 and the redesigned airfoil 
at M ^ = 0.60, Re = 20 million, and C/ = 0.40 



NACA 1412 

REDESIGN 

CONSTRAINTS 

a 

1.78 

0.87 


c m 

- 0.027 

-0.074 

-0.080 

c d 

0.0054 

0.0030 


x tr,u 

0.21 

0.59 

0.60 

x tr t l 

0.30 

0.48 (lam. sep.) 


tjnax 

0.120 

0.120 

0.120 

t at x = 0.20 

0.115 

0.096 

0.095 

t at x = 0.70 

0.073 

0.068 

0.070 

r ie 

0.0156 

0.0100 

0.0100 
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Table 2 contains a comparison of some characteristics of the redesigned airfoil with 
those of the NACA 1412 airfoil at the design flow conditions and lift coefficient The 
design constraints are also shown in this table for comparison. With the design method 
imposing the tolerances specified in Equations 68-72, the redesigned airfoil meets nearly 
all of the design constraints within the specified tolerances. Table 2 shows a 24 count 
reduction in drag due to the extent of NLF that was achieved on both surfaces. 

It took only six hours to complete the five iterations of the method that were required 
to redesign this airfoil. Approximately 20% of this time was associated with the Euler 
solver and CDISC airfoil design method, while 80% of this time was required by the sta- 
bility analysis code. 

To show that the final airfoil is nearly independent of the starting airfoil, the NASA 
High Speed NLF-0213 airfoil (ref. 8) was redesigned for the same flow conditions and 
design constraints as for the redesign of the NACA 1412 airfoil. 

The pressure distribution of the NASA High Speed NLF-0213 airfoil at = 0.60 , 

Re = 20 million, and c t = 0.40 is shown in Figure 27. The N-Factor envelopes for this 

pressure distribution and these flow conditions are shown in Figure 28. 

After 13 iterations of the NLF airfoil design method, the NASA High Speed NLF- 
0213 airfoil was successfully redesigned to meet nearly all of the imposed constraints. 
The analysis and target N-Factor distributions of the redesigned airfoil are shown in 
Figure 29, while the pressure distribution of the redesigned airfoil is shown in Figure 30. 
In Figure 31, the redesigned airfoil is compared with the NASA High Speed NLF-0213 
airfoil. 
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Figure 32 shows a comparison of the pressure distributions on the redesigned NASA 
High Speed NLF-0213 airfoil and on the redesigned NACA 1412 airfoil. Although they 
are not exactly the same, the pressure distributions do have similar shapes. The main rea- 
son that these pressure distributions are not identical is because the upper surface leading- 
edge pressures of the starting airfoils were very different. 

Table 3 was constructed to show the similarities between the characteristics of the two 
redesigned airfoils. With the exception of the angle of attack required at the design condi- 
tion and the pitching moment coefficient, both airfoils appear to have identical characteris- 
tics. If the pitching moment tolerance c m tol specified in Equation 69 had been reduced, 
perhaps the pitching moment coefficients of the two airfoils would be more similar. 


Table 3. A comparison of the design constraints and the 
characteristics of the NASA High Speed NLF-0213 and the two 
redesigned airfoils at = 0.60, Re = 20 million, and c/ = 0.40 



NLF-0213 

NLF-0213 

REDESIGN 

NACA 1412 
REDESIGN 

CONSTR. 

a 

1.03 

-0.18 

0.87 


c m 

-0.014 

-0.081 

-0.074 

-0.080 

Cd 

0.0042 

0.0030 

0.0030 


x tr,u 

0.30 

0.59 

0.59 

0.60 

x tr.l 

0.69 (lam. sep.) 

0.48 (lam. sep.) 

0.48 (lam. sep.) 



0.132 

0.120 

0.120 

0.120 

t at x = 0.20 

0.110 

0.096 

0.096 

0.095 

t at x = 0.70 

0.092 

0.068 

0.068 

0.070 

r le 

0.0095 

0.0100 

0.0100 

0.0100 
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In addition, the final target N-Factor envelopes are different, even though they both 
force transition to occur near 60% chord (see Figures 24 and 28). This may have also 
been due to the difference in the leading-edge pressures between the two starting airfoils. 

Figure 33 shows a comparison of the two redesigned airfoils. They appear to be very 
different. This is the case since the airfoil design method maintains the trailing-edge ordi- 
nates of the starting airfoil throughout the design process. As a result, to better compare 
the airfoils, the redesigned NASA High Speed NLF-0213 airfoil was rotated to match the 
average of the trailing-edge ordinates of the redesigned NACA 1412 airfoil. Figure 34 
compares the rotated NASA High Speed NLF-0213 redesigned airfoil with the redesigned 
NACA 1412 airfoil. In this figure, the similarities between these two airfoils can be more 
easily seen. 

In this plot, it appears that the two airfoils have a different leading-edge radius. The 
leading-edge radius is calculated by fitting a polynomial through the five points that com- 
prise the leading-edge, with these five points being ahead of 0.5% chord. In this region, 
the airfoils match very well, but then become different in the region between 2% and 20% 
chord. 

It should also be mentioned here that a change in a needs to be made to account for 
rotating the airfoil. To do this, a would have to be increased by 0.80°, which would 
increase the a of the redesigned NASA High Speed NLF-0213 airfoil at the design condi- 
tion from -0.18° to 0.62°. This better compares with the angle of attack at the design 
condition of the NACA 1412 redesigned airfoil, which was 0.87° . 
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5 .3 Airfoil for a Subsonic Transport Aircraft 

As a final example, the NASA Supercritical SC(2)-0412 airfoil (ref. 28) was rede- 
signed for the following flow conditions and constraints: 

M„ = 0.76 

Re = 10 million 
x lKU = 0.55 
Ci = 0.50 
c m = -0.100 

t max = 0.110 
t = 0.100 at x = 0.20 
t = 0.065 at x = 0.70 
r u = 0.0150 

A successful supercritical design is much more difficult to achieve than a subcritical 
design. With supercritical designs, if the specified target N-Factor distribution is not real- 
istic, then it is not possible to design an airfoil that has the desired target pressure distribu- 
tion, as was discussed in Section 3.1.2. As a result, trying to find a realistic target N- 
Factor distribution for a specific supercritical case is often a very tedious process. This 
was the case in this design. 

The pressure distribution of the NASA Supercritical SC(2)-0412 airfoil for the given 
flow conditions and design lift coefficient is shown in Figure 35. A shock is present on the 
upper surface near 40% chord. The upper and lower surface N-Factor distributions for 
this pressure distribution are shown in Figure 36. Using a transition N-Factor of eight for 
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the reasons discussed in Section 3.1.2, laminar separation occurred at the shock on the 
upper surface, while laminar separation occurred at 50% chord on the lower surface. 

Using the NLF airfoil design method, a new airfoil was designed in eight iterations. 
The upper and lower surface analysis N-Factors are plotted with the upper surface target 
N-Factors in Figure 37. Using the target N-Factor distribution shown in Figure 37, the 
pressure distribution shown in Figure 38 was calculated. Notice how the shock has been 
moved aft to 55% chord, and appears to be much weaker. Using the airfoil design method, 
the airfoil shown in this figure was then calculated. The redesigned airfoil is compared 
with the starting airfoil in Figure 39. 


Table 4. A comparison of the design constraints and the characteristics 
of the NASA Supercritical SC(2)-0412 airfoil and the redesigned airfoil 
at M m = 0.76, Re = 10 million, and C; = 0.50 



SC(2)-0412 

REDESIGN 

CONSTRAINTS 

a 

1.20 

0.58 


c m 

-0.075 

-0.101 

-0.100 

c d (wave) 

0.0017 

0.0002 


c d (total) 

0.0058 

0.0044 


x tr,u 

0.41 (shock) 

0.55 (shock) 

0.55 

x tr,l 

0.53 (lam. sep.) 

0.32 (lam. sep.) 


^max 

0.120 

0.107 

0.110 

t at x = 0.20 

0.109 

0.100 

0.100 

t at x = 0.70 

0.073 

0.065 

0.065 

r le 

0.0222 

0.0150 

0.0150 


A comparison of the characteristics of the NASA Supercritical SC(2)-0412 and the 


redesigned airfoil are shown in Table 4. With the exception of the maximum thickness 
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constraint, the final airfoil meets all of the design constraints imposed, even though the 
pitching moment constraint was released after five iterations. 

Since the laminar boundary layer separated at 32% chord on the lower surface of the 
redesigned airfoil, the viscous drag was not reduced in the design process. Nevertheless, 
due to the reduction in wave drag, the total drag of the redesigned airfoil was 14 counts 
less than that of the starting airfoil. Not only is this a result of moving the shock aft while 
maintaining the same lift coefficient, but it also results from allowing the N-Factors to 
grow as much as possible without forcing transition until the shock. 

The fact that the boundary layer remained attached on the lower surface only to 32% 
chord is a result of the geometric constraints that were imposed on the airfoil. As was 
mentioned previously, the lower surface target pressures are modified to meet the geomet- 
ric constraints. As a result, if the front spar thickness had been reduced to 9 or 9.5%, then 
perhaps the airfoil design method would not have had to work so hard to increase the pres- 
sures aft of 25% chord in order to meet the maximum thickness constraint. On the other 
hand, if the front spar thickness had not been reduced, then increasing the desired maxi- 
mum airfoil thickness to 12% would have given the same effect. 
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6.0 Concluding Remarks 

An automated two-dimensional method has been developed for designing NLF air- 
foils, while maintaining several other aerodynamic and geometric constraints. The 
method has been shown to work for a range of Mach numbers, Reynolds numbers, and air- 
foil thicknesses. The method has also been demonstrated for a supercritical case where a 
shock wave is present on the upper surface of the airfoil. 

In order to develop this NLF airfoil design method, several existing CFD codes were 
coupled together. In addition, a process was developed for calculating a target N-Factor 
distribution that forces transition to occur at the desired location. Using this target N-Fac- 
tor distribution, as well as the current analysis N-Factors and pressures, a method was also 
developed for calculating a target pressure distribution. Using this target pressure distri- 
bution, the current airfoil is redesigned to obtain a new airfoil that is closer to meeting the 
desired NLF, aerodynamic and geometric constraints. This method has been used to 
design a number of airfoils, with results shown for glider, commuter and subsonic trans- 
port applications. 

One advantage of this method is that it is capable of designing an airfoil in a short 
amount of time. Since an Euler solver has been coupled together with a turbulent bound- 
ary layer method to calculate the pressures over the airfoil, the design time is much less 
than that required for Navier- Stokes codes. As a result, a new airfoil with a large extent of 
upper surface NLF can be designed in only a few hours. 

In addition to the reduction in computer time required, a stability analysis code has 
been used to calculate N-Factors which are correlated to the transition location. Stability 
analysis methods have gained respect in the past few years and the prediction of the transi- 
tion location that results is taken as being fairly accurate. In this method, the stability 
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analysis code is automated to calculate the N-Factor distribution by varying the frequency 
of the T-S waves while assuming that the disturbances grow only in time. Since an N-Fac- 
tor distribution is calculated to determine the transition location of the airfoil, a design 
philosophy is presented for specifying a target N-Factor distribution for both subcritical 
and supercritical airfoil designs. Subcritical target N-Factor distributions are specified so 
that the flow does not undergo transition at slightly off-design conditions and reduces the 
uncertainty of the transition location by forcing the N-Factors to grow rapidly through the 
desired transition location. In supercritical designs, a target N-Factor distribution is spec- 
ified that forces the flow to transition before the shock so that laminar flow is not termi- 
nated at the shock. 

In order to design a new airfoil that possesses the desired target N-Factor distribution, 
an N-Factor-Target Pressure relationship was developed. This N-Factor design method 
relates a change in N-Factor at an x-location to a change in the local pressure gradient. In 
addition, this method is independent of Mach number and Reynolds number. 

Another attribute of the method is that it is capable of maintaining several aerody- 
namic and geometric constraints. A method was established to meet these constraints 
while also maintaining the desired amount of NLF on the upper surface of the airfoil. The 
approach implemented to meet these aerodynamic and geometric constraints is new. The 
method dictates that the upper surface target pressures are modified to meet only the NLF 
and aerodynamic constraints while the lower surface target pressures are modified to meet 
only the geometric constraints. 

This method has also been shown to be robust If enough design constraints are 
imposed, the airfoil that results is largely independent of the starting airfoil. Another 
advantage of this method is that the codes used have been coupled together in modular 
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form. This allows for other codes to be used in the place of any of the current compo- 
nents. The NLF airfoil design method works efficiently and well to design new NLF air- 
foils. Only by using this method could one appreciate how great it really works. 

There are several possibilities for extension of this research. The method could be 
applied to bodies other than airfoils and wings, with possible applications including fuse- 
lages and nacelles. In addition, the method could be extended to the design of airfoils for 
supersonic applications. Since large sweep angles are needed for supersonic wings, cross- 
flow instabilities would be a major issue. In these cases, boundary layer suction and blow- 
ing is often necessary to help reduce the crossflow disturbances. As a result, when 
extending the method to include supersonic designs, the method may also have to be mod- 
ified to account for suction and blowing. 
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Figure 1. A flowchart of the NLF airfoil design method 
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Figure 3. A flowchart of how COSAL is used to calculate the 
analysis N-Factor distribution for a pressure distribution 
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Obtain target N-Factor distribution 



Figure 4. A flowchart of the Target Pressure Design module 
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Figure 5. A typical subcritical target N-Factor distribution 
calculated by the Target Pressure Design module 



Figure 6. An example of an undesirable target N-Factor distribution 
for both subcritical and supercritical airfoil designs 
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Figure 9. The two regions that the upper surface pressures are 
divided into to calculate the target pressure distribution 




TARGET 
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(a) The first intermediate recovery pressures 



(b) The second intermediate recovery pressures 


Figure 11. The evolution of the target pressures in the recovery region 
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(c) The final target pressures in the recovery region 


Figure 11. Concluded 



Figure 12. The three regions that the upper surface target pressures are 
divided into in order to obtain the lift and pitching moment coefficients 
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Figure 13. A flowchart of the method by which the target 
pressures are modified to meet the aerodynamic constraints 
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Figure 14. The upper surface pressure distribution that is 
used to increase the leading-edge radius of the airfoil 



(a) Linearly scaling the upper surface target pressures 
Figure 15. The evolution of an initial lower surface target pressure distribution 
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(c) The initial target pressure distribution 


Figure 15. Concluded 
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Figure 16. A flowchart of the module for Modifying 
Target Pressures to Enforce Constraints 
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Figure 17. The pressure distribution and shape of the NACA 
64 r 212 airfoil at = 0.10, Re = 3 million, and c; = 0.30 



Figure 18. The upper and lower surface N-Factor distributions of the 
NACA 64 r 212 airfoil at = 0.10, Re = 3 million, and c, = 0.30 
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Figure 19. The upper and lower surface N-Factor distributions of the 
redesigned NACA 64p212 airfoil at = 0.10, Re = 3 million, and c/ = 0.30 
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Figure 20. The pressure distribution and shape of the redesigned 
NACA 64 r 212 airfoil at A/„ = 0.10, Re = 3 million, and C/ = 0.30 
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Figure 21. A comparison of the NACA 64j-212 airfoil 
and the redesigned NACA 64 t -212 airfoil 
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Figure 23. The upper and lower surface N-Factor distributions of the 
NACA 1412 airfoil at = 0.60, Re = 20 million, and c t = 0.40 
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Figure 27. The pressure distribution and shape of the NASA High 
Speed NLF-0213 airfoil at AC = 0.60, Re = 20 million, and c, = 0.40 



Figure 28. The upper and lower surface N-Factor distributions of the NASA 
High Speed NLF-0213 airfoil at AC = 0.60, Re = 20 million, and c, = 0.40 
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Figure 29. The upper and lower surface N-Factor distributions of the 
redesigned NLF-0213 airfoil at = 0.60, Re = 20 million, and c t = 0.40 
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Figure 30. The pressure distribution and shape of the redesigned 
NLF-0213 airfoil at = 0.60, Re = 20 million, and c; = 0.40 








Figure 31. A comparison of the NASA High Speed NLF-0213 
airfoil and the redesigned NLF-0213 airfoil 



Figure 32. A comparison of the pressure distributions of the 
redesigned NACA 1412 airfoil (Figure 25) and the redesigned 
NASA High Speed NLF-0213 airfoil (Figure 30) 
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Figure 33. A comparison of the redesigned NACA 1412 airfoil (Figure 25) 
and the redesigned NASA High Speed NLF-0213 airfoil (Figure 30) 



Figure 34. A comparison of the redesigned NACA 1412 airfoil and 
the redesigned NASA High Speed NLF-0213 airfoil after rotation 
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Figure 35. The pressure distribution and shape of the NASA Supercritical 
SC(2)-0412 airfoil at AC = 0.76, Re = 10 million, and c, = 0.50 



Figure 36. The upper and lower surface N-Factor distributions of the NASA 
Supercritical SC(2)-0412 airfoil at AC = 0.76, Re = 10 million, and C/ = 0.50 
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Figure 37. The upper and lower surface N-Factor distributions of the 
redesigned SC(2)-0412 airfoil at = 0.76, Re — 10 million, and c ; = 0.50 



Figure 38. The pressure distribution and shape of the redesigned 
SC(2)-0412 airfoil at = 0.76, Re = 10 million, and C[ = 0.50 
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Figure 39. A comparison of the NASA Supercritical SC(2)-0412 airfoil 
and the redesigned SC(2)-0412 airfoil 
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Appendix A. A Linear Scaling Method 


In a few cases within the NLF airfoil design method, it is desirable to calculate a new 
target pressure distribution that has nearly the same shape as another pressure distribution, 
but different pressure magnitudes. This can be accomplished by linearly scaling an exist- 
ing pressure distribution. Consider the pressure distribution shown in Figure A.l. Sup- 
pose that this pressure distribution, C p . , is to be linearly scaled to obtain a new target 

pressure distribution. In addition, assume that the pressure coefficient at 50% chord of the 
new target pressures is to be -0.30, while the leading-edge pressure coefficient is to remain 
unchanged. The new target pressures can be calculated using the relation 

c p.T.j= <- c p.r c P ,0 G+c r.i (A*) 

where G is the scale factor and C p j represents the leading-edge pressure coefficient of 
C . . In this case, 

p.j 


- 0.30 - C„ 

G = p — 

C -C 

S. * S. i 


(A.2) 


where k represents the station nearest to 50% chord. Notice that when Equations A. 1 and 


A.2 are used, C p T j = C p j and C p T k = -0.30 as desired. The new target pressure 
distribution that results from using Equations A. 1 and A.2 is shown in Figure A.2. 
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Figure A.l. The pressure 
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distribution used to calculate the new target pressures 



Figure A.2. The new target pressure distribution 
that results from using a linear scaling technique 
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Appendix B. A Weighted Averaging Technique 


On several instances within the NLF airfoil design method, a weighted averaging tech- 
nique is used to calculate a new target pressure distribution from two existing pressure dis- 
tributions. This weighted averaging technique calculates a new target pressure 
distribution using the relation 


C P,TJ= W j C p.j,l + V- W ? C PJ.2 (B-D 

where W. is a weighting function, and C p ■ x and C p ■ 2 are two existing pressure distri- 
butions. A similar weighted averaging technique is used in reference 29. 

In the NLF airfoil design method, Equation B.l is most useful when Wj is allowed to 

vary along the chord since it is usually desirable to maintain certain characteristics of each 
of the two existing pressure distributions. Consider the two pressure distributions shown 
in Figure B.l. Suppose that these two distributions are to be used to calculate a new target 
pressure distribution, and that the leading-edge pressures of C p . , and the trailing-edge 

pressures of C . 2 are to be retained as the target pressures in each of the respective 

regions. If W- is allowed to vary from a value of 1 at the leading-edge to a value of 0 at 

the trailing-edge, it can be seen from Equation B. 1 that the resulting target pressures 
would have the desired properties in the leading-edge and trailing-edge regions. 

The simplest expression for W. that satisfies these requirements would be a linear 
variation between the leading-edge and the trailing-edge. That is. 




(B.2) 
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Figure B.2 shows the new target pressures that result from using C p j j and C p - 2 from 

Figure B.l and the weighting function from Equation B.2. 

Equation B. 1 can also be used to modify the taiget pressures over only a small region 
of the chord. In fact, this is the only application of the technique that is used within the 
NLF airfoil design method. Suppose that C p j 1 represents the current target pressures 

and C p j 2 represents the current analysis pressures. It may be desirable to retain the tar- 
get pressures of C p j j ahead of 60% chord (i.e., Wj = 1 ahead of 60% chord), but, at the 
same time, retain the characteristics of C„ , , at the trailing-edge. If station k represents 

r* J* * 


the ordinate closest to 60% chord, then Wj can be varied linearly from a value of 1 at sta- 
tion it to a value of 0 at the trailing-edge. As a result, the weighting function aft of 60% 
chord can be represented by the expression 



iis 

1 -X k 


(B.3) 


The target pressure distribution that results from using Equation B. 1 is shown in Figure 


B.3. 
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Figure B.2. The new target pressure distribution obtained by 
using a weighting function over the entire length of the chord 
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Figure B.3. The new target pressure distribution obtained by 
using a weighting function over a small region of the chord 
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